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CALCULATION OF CAMBER USING SLENDER WING THEORY 
AT MACH NUMBER 1.0* 


By William B. Igoe 
Langley Research Center 

SUMMARY 

A calculation procedure based on slender wing theory at Mach number 1.0 is pre- 
sented for the determination of the mean line camber ordinates for laterally symmetrical 
wings with polygonal planforms. A uniform chordwise load and a continuous but otherwise 
arbitrary spanwise load distribution is assumed. The procedure is generally limited to 
slender wings at high subsonic speeds and low lift coefficients. 

Graphs of functions required in the calculation procedure are presented. In addition, 
an application of the method of calculation for a variable- sweep wing planform is described 
and some calculated camber results are shown. Although the slender wing theory for 
Mach number 1.0 is considered to be inadequate for the calculation of the detailed shape of 
airfoil section camber for supercritical wings designed for near-sonic speeds, a compari- 
son with experimental results shows that it predicts the general magnitude and trend of the 
spanwise variation of the wing twist angle. 

INTRODUCTION 

Slender wing theory was introduced in 1946 by Robert T. Jones in reference 1, and 
in 1948, the application to sonic flow was stated explicitly by Heaslet, Lomax, and Spreiter 
in reference 2. At about the same time, similar results were obtained in England by 
A. Robinson and in France by R. Legendre, E. A. Eichelbrenner, and A. von Baranoff. 
Detailed references to these developments are contained in volumes 6 and 7 of the 
Princeton Series on High Speed Aerodynamics and Jet Propulsion. (See refs. 3 and 4.) 

As shown in reference 2, the linearized potential flow theory equation relating the 
downwash of a slender lifting wing to its load has an especially simple form at Mach 1. 

This equation for the downwash at a field point on a wing at Mach 1 has the same form as 
the incompressible flow equation for the downwash in the far wake (neglecting rollup of the 
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trailing- vortex sheet) downstream of the field point for that part of the wing ahead of the 
field point. The relationship between the sonic and incompressible flow fields is easily 
visualized when the Prandtl-Glauert compressibility factor for Mach 1 is considered as 
stretching the wing dimensions infinitely in the stream wise direction. The pressure sig- 
nals in the zone of silence behind a sonic Mach line through the field point cannot propa- 
gate upstream in sonic flow. Therefore, only that part of the wing forward of the field 
point contributes to the downwash at the field point. 

The simple form of the downwash equation for a slender lifting wing in sonic flow and 
its resemblance to the far-wake downwash equation of incompressible flow allows a simple 
solution to be used. This solution is similar to the lifting- line theory solution of Glauert 
in reference 5 for incompressible flow. By following Glauert's method, the sonic wing 
downwash velocity, and thereby the local stream wise slope of the wing- section cambered 
mean lines, may be expressed in terms of the coefficients of a Fourier sine series. This 
series represents the partial spanwise load distribution of that part of the wing ahead of 
the field point. The local cambered mean line slopes may then be integrated stream wise 
to obtain cambered mean line ordinates. 

For the special case of a polygonal wing planform with lateral symmetry and with 
uniform chordwise load, the Fourier coefficients for computing cambered mean line slopes 
and similar coefficients for computing cambered mean line ordinates have been derived 
for continuous but otherwise arbitrary spanwise load distributions. The sum of the 
resulting series expression for cambered mean line ordinates has been taken to obtain 
a closed-form solution. This closed-form solution avoids the necessity of evaluating 
a term-by-term summation of a slowly convergent series which requires a large num- 
ber of terms for numerical accuracy. The resulting solution is applicable to a variety 
of wing planforms, including those typical of variable- sweep wings. It can be extended 
to other planforms by means of superposition. Ths usefulness of the slender wing 
theory at Mach 1 is mainly for wings at high subsonic speeds. As indicated in refer- • 
ence 2, this simplified form of the linear theory should be applicable for wings for which 
the slenderness parameter j3A is small. 

Graphs of the functions encountered in the closed-form solution are presented and 
an application of the method of calculation for a variable-sweep wing planform is 
described. Some calculated camber results are included and the theoretical wing twist 
distribution for a supercritical wing designed for near-sonic speeds is compared with an 
experimentally obtained twist distribution. 
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SYMBOLS 


A aspect ratio 

A m coefficient in trigonometric sine series representing the wing- section 

design lift coefficient distribution (see eq. (A4)) 

B n coefficient in Fourier sine series representing the partial wing span wise 

load distribution (see eq. (4)) 

b wing span 

Cl wing lift coefficient 

^n coefficient in series solution for wing mean line camber ordinates 

(see eq. (8)) 

Cp pressure coefficient 

^P (pp) lower surface (^p) upper surface 

c wing chord (stream wise, in general, except for ct with a < 

c av average chord, S/b 

c i section lift coefficient 

D m function in analytical summation of series (see eq. (Cl)) 

dj perpendicular distance from wing pivot point to wing leading edge 

d£ perpendicular distance from wing pivot point to wing trailing edge 

d 3 perpendicular distance from wing pivot point to wing -tip chord 

E m function in analytical summation of series (see eq. (C4)) 

F,G,H,I functions of n, m, and 6 (see eqs. (All), (A12), (B6), and (B7)) 
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K m function in recursion formula solution for A m coefficients for elliptical 

spanwise load distribution on straight- tap er ed wings (see eq, (E5)) 

k = tan A (with exception of k 3 = -tan a) 

k* = tan A* (with exception of kg = -tan a') 

L m,s functions in summation for wing mean line camber ordinate (see eq. (14)) 

l length of variable- sweep wing in fully swept position with straight leading 

edge (before truncation) 

M free-stream Mach number 

m index in trigonometric series for section lift coefficient distribution 

(odd integers only for wings with lateral symmetry) 

n index in Fourier series for P* (odd integers only for wings with 

lateral symmetry) 

P* partial wing spanwise load distribution function (see eq. (2)) 

p integer 

r index of summation 

S wing area 

S* partial wing area 

s index of summation of L mjS functions (see eq. (14)) 

V free-stream velocity (parallel to X-axis) 

w downwash velocity (parallel to Z-axis) 

x,y,z stream wise, spanwise, and vertical orthogonal Cartesian coordinates with 

origin at the wing leading- edge apex 
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streamwise distance from leading- edge apex to intersection with X-axis 
of line projected from outboard leading edge of variable-sweep wing panel 


Xb 


x le 


streamwise distance from leading- edge apex to intersection with X-axis of 
line projected from wing tip chord jfor a < -| j 

streamwise location of leading edge of local chord 


Axq local streamwise distance from wing leading edge to field point 

Yl span of variable-sweep wing in fully swept position (before truncation) 

yj spanwise point at which Mach line (xq = Constant) through field point 

intersects the wing leading edge 


Y 2 spanwise point at which Mach line (xq = Constant) through field point 

intersects the wing trailing edge 

yg spanwise point at which the Mach line (xq = Constant) through field point 

intersects the inward slanted wing tip (tor a < 

Azq camber ordinate measured with respect to leading edge of local streamwise 

chord : (Az 0 = z 0 - z le ) 

j3 = (1 - M 2 ) 1 / 2 


A difference 


8 angle between wing outboard leading and trailing edges (d = Aj - Ag'j 

e wing twist angle 


6 = cos~l(-y) 

A sweep angle 

X taper ratio for straight- tapered wings 

a angle of wing tip chord with respect to Mach line fa = o 4 + A| - Aj 
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Subscripts: 


A,B lower and upper limits, x A corresponds to wing leading edge and Xg 

corresponds to wing trailing edge or line xq in general, y A and yg 
correspond to left and right wing tips or to -yj . and yj, respectively. 

# A and are similarly defined 

k wing leading- edge kink for variable- sweep wings 

le leading edge 

p wing pivot 

r wing root 

te trailing edge 

t wing tip 

0 field point 

1 wing leading edge 

2 wing trailing edge 

3 wing tip 

Superior or overhead bar ( — ) indicates dimensional length. Unbarred lengths are 
nondimensionalized with respect to the wing semispan. Asterisk (*) superscript indicates 
partial wing parameter. Prime (’) superscript indicates wing geometric characteristics 
when leading edge is straight for variable-sweep wing. 

BASIC EQUATIONS 


In linearized wing theory, the equation relating the downwash or normal induced 
velocity at a field point to the partial loading ahead of the field point for a slender lifting 
wing at M = 1 is given by equation (62) of reference 2. This equation may be written as 


3P* 

wo 1 f^B 9y 
T " "477 Jy A y 0 - y 


(i) 
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where 


P* = f XB AC dx 
Jx A P 


( 2 ) 


and the limits of integration are established by the boundaries of the region S* as shown 
in figure 1. Satisfaction of the Kutta condition at the wing trailing edge is assured by 
imposing the condition that the load distribution is zero behind the trailing edge. The 
velocity ratio in equation (1) is equal to twice the normal induced velocity ratio in the 
near wake for incompressible flow (see ref. 5, for example) for a lifting- line wing con- 
sisting only of the region S*. The factor of two accounts for the difference between near- 
wake and far-wake conditions if rollup of the wing trailing vortex sheet is neglected. 

A well-known method of evaluating the improper integral of equation (1) with the use 
of Fourier analysis was shown by Glauert (ref. 5). By following reference 5, the spanwise 
variable 6 is introduced so that 

y = -cos 6 (3) 


where 6 varies from 0 to u as y varies from -1 to 1 as shown in figure 2. 
The partial load function P* is then expanded in a Fourier series 


p* 


■l 

n=l 


B n sin nd 


(4) 


With the use of equations (3) and (4), the integral of equation (1) becomes 


^0 

V 


j_ r e B 

4l Jf) A 


CO 


2 


nB n cos n0 


n=l 

cos 5" - cos 


d6 


(5) 


The order of summation and differentiation has been interchanged in obtaining equation (5). 
Reference 6 shows that this interchange of order is permissible if the function P* is 
continuous and if its derivative is sectionally continuous in the interval ~ir Ik 6 = it. In 
general, this requirement is satisfied for wings with span load distributions which are 
everywhere continuous. The conditions on P* are also sufficient to permit the inter- 
change of order of summation and integration. Accordingly, the integral of equation (5) 
may be taken inside the summation sign and evaluated as shown in reference 5 to obtain 

OO 

W Q _ 1 V nB n Sin n6 0 

T " " 4 L sin % W 

n=l 
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In this step the lower and upper limits of integration for 6 have been taken as 0 and tt, 
respectively. The extension of the limits for 8 beyond those defined by S* is possible 
if the function P* is defined as being zero outside the region S*. 

The slope of the mean line camber is equal to the down wash velocity ratio 


dz _ w 
dx V 

Xnte 


( 7 ) 


where 


Hence the slope of the mean line camber may be obtained directly from equation (6). 
gration of equation (6) with respect to x gives the mean line camber ordinate 


z = 


1 

4 sin 8q 


^ nC n sin n#Q 


n=l 


C 


n _ 



B n dx 


( 8 ) 


and again the order of summation and integration has been interchanged. By taking the 
integral limits on x from 0 to the field point xq, the camber ordinates are obtained 
with respect to the leading- edge apex. An indeterminate form is encountered in equa- 
tions (6) and (7) for field points at the wing tip where 8q = n. The limiting form is seen 

to be 


lim 

8 ^ ji sin 0 


n(-l) n+ ^ 


Thus, for $q = tt, equations (6) and (7) may be rewritten as 

OO 

^r=-f £ n2B n (-l) n+1 ( 6 0 = 7r) (9) 

n=l 


and 

OO 

z = - 1 £ n2 G n (-i) n+1 (e 0 = 7f) (10) 

n=l 

The series summations of equations (9) and (10) are generally not convergent so that 

w Q 

— — 00 and z — 00 as 6q — ir. Equations (6) and (7) have not been restricted as to wing 
planform or load distribution with the exception that the partial load function P* must be 
continuous and have a sectionally continuous first derivative. In the following analysis, 
the applications are restricted to wings with lateral symmetry and with a uniform 
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chordwise load distribution. This chordwise load distribution causes infinite camber 
slopes to appear at the leading and trailing edges but the camber ordinates are finite 
except at the discontinuities represented by the wing root, wing tip, and wing leading- edge 
sweep discontinuity in the case of a variable-sweep wing. 

For a uniform chordwise load distribution, the B n Fourier coefficients of equa- 
tion (4) are evaluated in appendix A for a wing planform and a field-point location similar 
to that shown in figure 2. For the same conditions, an evaluation of the C n coefficients 
of equation (8) is shown in appendix B. The results of an analytical summation of the 
series of equation (7) are presented in appendix C and an example of the procedure used 
in obtaining the summation terms is shown in appendix D. Some methods of obtaining 
the A m coefficients which are needed to define the spanwise section lift coefficient dis- 
tribution 

n = l sin rad (11) 

m 

are discussed in appendix E. 

The application of the method of calculating camber is shown for wings with less 
restricted planform geometry and field-point location than are considered in appendixes A 
and B. However, the limitation of uniform chordwise load distribution is retained. 

APPLICATION TO WINGS OF SPECIFIED PLANFORM 

The equations developed in the preceding section and in appendixes A, B, and C are 
extended and applied to a wing planform with discontinuous slope (sweep angle) of the 
leading edge. This planform configuration is encountered in variable-sweep wings of the 
type shown in figure 3. The leading edges are idealized as straight continuous line seg- 
ments. The trailing edges are taken as continuous straight lines with continuous sweep 
angles. Wings with swept forward leading or trailing edges or with wing tips that are 
leading edges are not considered. The leading- edge sweep of the outboard portion of the 
wing is limited to angles equal to or less than the leading- edge sweep of the fixed inboard 
part of the wing. The latter conditions may be expressed as 

0 ^ A 1 s A'j 

0 S A 2 ^ Aj - 6 

° = a = i 
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Although other planform configurations are possible in variable- sweep wings, the config- 
urations included under these conditions are considered the ones most likely to be encoun- 
tered in practice. 


Wing Geometry 

The equations in the preceding section and in appendixes A, B, and C have been 
developed for a wing of unit semispan. However, for variable- sweep wings the geometry, 
including the semispan, changes with sweep angle. It will be convenient in the camber cal- 
culations to consider certain geometrical parameters describing the variable-sweep plan- 
forms to be dimensional in character. The wing geometry will then be nondimensionalized 
in terms of the semispan for each wing- sweep angle. The dimensional parameters will be 
indicated by a superior or overhead bar. Plain or unbarred parameters will be nondimen- 
sional. The following fixed constants will be considered given: Aj, 6, a\ Xp, yp, ds, 
c t , and cr. In addition, the variable - sweep angle Aj of the leading edge of the outboard 
part of the wing will also be given. 

The quantities which are needed to introduce the effects of variable-wing geometry 
into the wing camber calculations can be computed from the input constants and Aj. The 
physical significance of most of the following geometrical quantities can be inferred from 
the sketch of figure 3. The first step is to obtain the angles and slopes of the various 
edges of the wing under variable- sweep conditions: 

A 2 = A 1 - 6 

A 2 = A i- 6 

a = a 1 + Aj - Aj 

k| = tan Aj 

f t 

kj = tan Aj 

k 2 = tan A 2 

i « 

kg = tan Ag 

k 3 = -tan a 

t . 

krj = -tan a 
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The next step is to obtain the length parameters which are needed to establish the nondi- 
mensional geometry and the criteria for determining field-point regions under variable- 
sweep conditions : 

d x = (x p - kjy p ^cos Aj 

b _ — 1 / ^3 \ 

2 y P + k! - kg \cos a + cos Ajy 

| = y p + S3 

El = y + 1 ( d 3 , d l \ 

2 P kj - kg y cos a ' cos A ij 

^ = vdrxr kl5?p 

_ _ d 3 _ 

x b = x P + ^i^ +y p tana 

d 2 = fc' r + kgy p - 5^ cos Ag 
_ _ d 2 

Cr = x p + ^rr 2 - k2y p 

x a 

y k = -7— 

kj - ki 

y k = cos 2 A t 1 (tjip + y^j 

*k = k iyk 

? ? ? 1 — — — ? . 

The computed quantities Ag, kp kg, kg, dj, d 2 , b/2 do not change with the 
variable-sweep angle. The remainder of the quantities are functions of the sweep angle. 
The quantities x a , xb, c r , y k , and x k are divided by the semispan b/2 to obtain 
the nondimensionalized values of these quantities. The dimensional semispan b/2 
divided by itself becomes the nondimensional semispan y^. = 1.0. 

The method of calculation is, of course, applicable to other than variable-sweep wing 
planforms. For instance, fixed s trai ght- tap e r ed sweptback wings can be represented by 
specifying a pivot axis which coincides with the wing leading- edge apex. Fixed wings with 


( A i < A i) 

( A 1 - A',) 
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discontinuous leading- edge sweep can be represented by specifying a pivot axis which coin- 
cides with the wing leading- edge sweep discontinuity. In the latter case, may not be 
known and can be determined by the following relation: 


Cr — 


cos h-2 
cos Ag 


c r + 


1 - 


cos Ag 
cos a: 


Pp - k 2 - k 2 


cos Ag 
cos A 2 


KP 


( 12 ) 


The wing geometry relations presented in this section are sufficient for a complete 
description of a variable-sweep wing as required for a calculation of the camber. 


Field-Point Locations 

In specifying field-point locations, it is convenient to use streamwise coordinates 
which are measured from the leading edge of the local chord and are expressed in frac- 
tions of that chord 

Axq _ xp - x le 
c c 

Here x^ e is the streamwise location of the leading edge of the local chord, measured 
from the wing leading-edge apex. For variable -sweep wings of the type described in the 
preceding sections, the field-point coordinates may be calculated with the aid of equa- 
tion (13) as follows : 

For 0 < yQ § yjj and 

A 

c = c r + (k 2 - kj)y 0 

x o = (nr) [ c r + ( k 2 - k i) yo] + k iyo 

For y^ < yo = Yt " c t cos a a ' = a = f j 

c = c r - xa + (k 2 - kj)y 0 
x 0 = ^-^jc r - xa + ^k 2 - kj) y 0 ] + xa + k^Q 
For y t - ct cos a < yo = yt and a’ £ a ^ 
c = x b - Xa + (k 3 - kj) y 0 
x () = (~?)[ x b " % + ( k 3 - k l) Yo] + % + k lYQ 
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The foregoing relations are not necessary for a satisfactory specification of field 
points. Any suitable method of specification may, of course, be used. One other method 
which was found useful was to specify field-point - coor dinate locations on the wing in the 
fully sweptback position (with the leading edge straight). At the forward sweep positions, 
the field points were considered to have moved with the wing. A description of some cam- 
ber results obtained in this way is presented in a subsequent section. 


Calculation of Mean Line Camber Ordinates 

In addition to the field-point coordinates (xg,y 0 ) and the kink coordinates (x^y^, 
the following quantities are used in the calculation: 


XQ 

71 "if 


Yl = 

= 

= 


X Q - 

ki 

x 0 ~ c r 
k 2 

X Q - x b 
k 3 


(° = x 0 S x k) 
(x k Sx 0 5x a + k x y t ) 

(c r = xq] 

x a + kjy t < xg and a < ■ 


From these values of y, the values of 6 are obtained from equation (3), in general, as 


6 = cos _1 (-y) 

so that, for example, 9q = cos~-^-y^ and % = cos"l^-yj £ ^, etc. 

By rewriting equation (7) with the use of equations like equation (B5) for C n and 
substituting the summation equations (Cl) to (C4), the expression for the mean line cam- 
ber ordinate becomes 


z 0 


7r sin 9q 


A 


m 


m,s 


m 


( 14 ) 


The summation on m is made over all values of m for which the A m coefficients are 
known. The L m g functions are as follows: 
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- k| - k 2 |lS m (0i) - E m ^ 

+ |(k| cos 2 0 2 - k| cos 2 6 1 )D m (7r) - (k| cos 0 2 " kf cos (15) 

L m,2 = 5^0 - (*a + k iyt J 2 D m (7r) + k£x 0 - (x* + k iy JJ ^D^tt) - D m (|) (16) 

■^m,3 = "k 2 E m(^3) + k 2 E m ^j “ (k| “ k^EinOr) 

~ 4^1 " k 2 sin 2 %)D m (ir) - jkf - k|(l + cos ®3)]Dm(^ (17) 

L m,4 = (k'l - ki)(x 0 - x k ) D m (%) - D m ^j + |(cos %)D sn (u)j 

- (k'j 2 - k 2 )£ m (%)- E m (|) + I(cos 2 0 k )D m OO - (cos %) (18) 

L m,5 = k| E m (^ 2 ) - E m ^ - (k| cos &2j D m(^J “ ^ os @2) D m( 7r ) (19) 

The particular L m g functions to be used and the values of 6 to be inserted in these 
functions are specified for the various ranges of the field-point locations. 

(1) For field points forward of the wing leading- edge kink and the trailing- edge 
root (0 < XQ £ x k ^ c r or 0 < xq = c r ^ x^J : 

l"m,s = km,l 

s 
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with 


Q\ = cos 





where kj is used in L m ^ in place of k]_. 

(2) For field points forward of the wing leading- edge kink but behind the wing 
trailing- edge root fc T < xq = x k \: 


l 


^m,s _ ^m,l 


with 


01 = cos 




0=1 = cos - * 


c„ - 


X 0 


k 2 


and again kj is used in L mj i in place of k l- 

(3) For field points behind the wing leading- edge kink but forward of the trailing- edge 
root and of the outboard leading- edge tip ^x k < xq = c r = % + k^y^ or 
x k < x Q S xa + k lYt g c r ): 


I 


km,s _ k m,l + ,4 
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with 


6\ = cos 


-1 *a I x 0 

k l 


£ 

2 


% = cos_ 1 (-yk) 

(4) For field points behind the wing trailing- edge root and of the wing leading- edge 
kink but forward of the wing leading- edge tip ^ c r < xg = x a + kjy t or 

c r= x k <x 0- x a + k iyt) 

L m,s = L m,l + L m,4 


„ _1 x a - x 0 

01 = COS 1 

1 k l 


09 = COS 


-1 c r " x 0 


k 2 


a 


7T\ 


or 


% - cos _1 (-y k ) 

(5) For field points behind the wing leading- edge tip and the trailing- edge root, with 
j streamwise wing tips 
% + k 'iy t = c r < x 0 = c r + k 2yt) : 

L m,s = L m,l + L m,2 + L m,4 


treamwise wing tips (c T = x a + k^y t < Xq S c r + k 2 y^ 


with 


= 71 


02 = cos: 


-X c r ~ x 0 


0 k = cos-^-yk) 
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( 6 ) Same as (5) except fct r S a < ~) wing tips not streamwise 

x h - c r \ z/ x b - c \ 

. i. „ ^ c D r „ . i,_*, < „ _y L\ . 


c r =xa + kjyt < x 0 S 


k2 - ks 


or ^ + kiy t S c r < x 0 S 


k2 - k 3 J 




L m,2 + L m,3 + L m,4 + L m,5 


i L r 

6*2 = cos A — 




k 2 

J X b - XQ 


= cos ^ 


% = cos _1 (-yk) 


The camber equations given in this section complete the presentation of the method 
of calculating camber for variable-sweep wings using the sonic slender wing theory. The 
wing geometry has been specified, field-point locations have been determined, and the cam- 
ber relations for various regions of the wing have been stated. The way in which the 
method can be applied to wings of simpler planform has also been indicated. The next 
section will present camber results which have been obtained for some specific planform s. 


Calculated Camber Results 

Camber ordinates have been calculated for the variable- sweep planform shown in 
figures 3 and 4. At the design sweep angle of 75° 58’, the leading edge is straight and 
the wing has an aspect ratio of 1.64. At the most forward sweep of 15° 58% the wing tips 
are streamwise. Camber results for these two wing sweeps and for a third intermediate 
sweep angle of 45° 58® are presented in figure 5. 

The wing loading imposed for all sweep positions was elliptical spanwise except in 
the vicinity of the wing tips which were intentionally unloaded. Consequently, the ordi- 
nates presented for each wing- sweep position represent the camber required to achieve 
nearly elliptical span load distribution at that wing sweep. As explained in a preceding 
section, the variable- sweep wing planform has been idealized in the region of the pivot by 
extending the inboard and outboard leading edges to a common intersection at the kink 
point and by extending a straight trailing edge into the center line for all wing- sweep 
positions. 

The camber ordinates are calculated with respect to an orthogonal Cartesian coordi 
nate system with its origin at the leading- edge root of the fixed inboard part of the wing. 
These ordinates are therefore absolute in that no shear or dihedral variation has been 
applied to them. The field points (points at which the camber has been calculated) were 
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located along constant x lines for the wing in the fully swept (75° 58') position. For the 
other two sweep positions, the field points were maintained in the same relative positions 
with respect to the local planform geometry so that those points which were located on the 
outer wing panel moved with the panel as the sweep was varied. (See fig. 4). The cam- 
ber ordinates in figure 5 have been plotted as a function of the original field-point positions 
on the variable-sweep wing in the fully swept position. By this method the calculated cam- 
ber lines for the wing can be traced through the various sweep positions. Of course, only 
a deformable wing could satisfy the conflicting camber requirements for the various sweep 
positions. However, by a suitable choice of pivot axis orientation and initial wing shear or 
dihedral variation in the fully swept design condition, the conflicts may be minimized and 
a satisfactory compromise may be possible. This aspect of the problem has not been 
treated herein. For comparison purposes, supersonic (M = 3) camber lines are also 
shown in figure 5 for the 75° 58’ sweep angle. The supersonic camber was calculated by 
using references 7 to 9. 

Mean line camber ordinates were also calculated for a conventional sweptback wing 
with straight taper and stream wise tips. The wing had an aspect ratio of 8, a taper ratio 
of 0.3, and a quarter- chord- line sweepback angle of 40°. As mentioned previously, wings 
of this type can be represented in the equations of the preceding sections by specifying a 
pivot axis which coincides with the wing leading- edge apex. Camber ordinates are shown 
in figure 6 for span wise locations at 0.1, 0.3, 0.5, and 0.8 of the semispan and are com- 
pared with ordinates (presented in ref. 10) which were calculated by a different method. 

In figure 6, the camber ordinates have been adjusted to give a zero camber ordinate along 
the leading edge. The wing section of the model of reference 10 was formed by distributing 
standard symmetrical airfoil thickness profiles (NACA 65A-series sections) from refer- 
ence 11 about the calculated mean line camber ordinates. 

A limited appraisal of the value of slender wing theory for design purposes may be 
obtained from the aerodynamic data of reference 10. The data indicate that some modifi- 
cations to the basic wing design were necessary to improve aerodynamic performance in 
the high subsonic speed range. The primary modification was an extension of the inboard 
leading edge. This extension increased the sweep angle and reduced the thickness ratio 
and leading- edge camber in that region. The local supersonic flow on the inboard portion 
of the wing was relieved and the formation of strong shocks which were encountered in that 
region was delayed by the modification. Further modifications to wings of this type were 
reported in reference 12 and confirmed the results shown in reference 10. 

The wing planform of reference 12 resembled the planform of variable- sweep wings 
in an intermediate- sweep position, but it was a fixed- geometry wing configuration. As 
previously mentioned, wings of this type can be represented in the equations of the 
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preceding section by specifying that the pivot axis coincides with the wing leading- edge 
kink position and by employing equation (12) to relate c^ to c r . 

The slender wing (linear) theory for Mach 1 is not considered to be adequate for 
the calculation of the actual detailed shape of the airfoil section camber for supercritical 
wings designed for near-sonic speeds. These wings tend to be heavily loaded and have a 
large region of supercritical flow. Consequently a more elaborate nonlinear transonic 
analysis in which a separation of camber and thickness effects is not possible or a trial- 
and- error experimental development is considered to be necessary. Despite this limita- 
tion, an application of the linear theory for the prediction of wing twist distribution for 
wings of this type has been examined. 

The results of a twist calculation for a supercritical wing planform designed for 
near-sonic speeds is presented in figure 7 where the twist angle e was calculated from 



The experimental data shown in this figure are based on unpublished results obtained at 
the NASA Langley Research Center by Charles D. Harris, Charles H. Fox, Jr., Robert T. 
Taylor, Dennis W. Bartlett, and Richard J. Re. An elliptical spanwise load distribution 
was assumed for the calculation and was approximately the same as the measured span- 
wise load distribution. The experimental twist distribution was obtained from structural 
wing twist data and was corrected for aeroelastic effects. An adjustment for wing inci- 
dence and the effect of angle of attack was also made. 

A comparison of the twist distribution predicted by the linear theory with that 
obtained experimentally shows that the linear theory predicts the general magnitude and 
trend of the spanwise variation of the twist angle. However, for a unit lift coefficient 
(Cl = 1.0), the experimentally determined twist distribution exceeds the calculated results 
by about 2° to 4° (excluding the immediate region of the kink). An experimentally deter- 
mined adjustment of the wing incidence angle may therefore be necessary in applying the 
theory. It should be noted that the experimental results include the effect of a fuselage 
extending approximately over the inboard 10 percent of the wing; however, fuselage effects 
were not included in the linear theory prediction. 

CONCLUDING REMARKS 

A calculation procedure based on slender wing theory at Mach 1 has been derived for 
the determination of the mean line camber ordinates for polygonal planform wings with 
lateral symmetry. A uniform chordwise load and a continuous but otherwise arbitrary 
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spanwise load distribution was assumed. The range of applicability of the results is gen- 
erally restricted to low lift coefficients and high subsonic speeds for slender wings for 
which the slenderness parameter 0A is small. 

Although the slender wing theory for Mach 1 is considered to be inadequate for the 
calculation of the detailed shape of airfoil section camber for supercritical wings designed 
for near-sonic speeds, a comparison with experimental results shows that it predicts the 
general magnitude and trend of the spanwise variation of the twist angle. An experimen- 
tally determined adjustment of the wing incidence angle may, however, be necessary in 
applying the slender wing theory to supercritical wings. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 28, 1971. 



APPENDIX A 


EVALUATION OF THE B n FOURIER COEFFICIENTS 


The B n Fourier coefficients for a wing of arbitrary planform with a spanwise load 
distribution which is everywhere continuous are given by 

E n = | P* sin n6 dfi (Al) 


where P* has been assumed to be an odd function of 6. If the wing has lateral symme- 
try, that is, symmetry about the x,z plane at y = 0, then only odd harmonics will occur 
in the B n coefficients. Consequently, equation (Al) may be changed to an integral over 
the right half of the wing 




P* sin n8 d9 


(A2) 


'7J 

For this case P* is not only an odd function of 6 but is also symmetrical about 6 = 
The further restriction of uniform chordwise load distribution is now introduced so that 


^^p c l 

Consequently, equation (2) may be rewritten as 


P* = 



(A3) 


The spanwise distribution of c i may be expressed in the trigonometric series of 
equation (11) 


Ci = } A m sin m6 (A4) 

m 

The series for c i can be a finite trigonometric series or it can be a Fourier series. In 
the latter case, generally just the first few terms are needed for c i distributions which 
are smooth. The use of equation (A4) in equation (A3) yields 

P* = (x B - ^ A m sin m 6 (A5) 

m 
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APPENDIX A — Continued 


Equation (A2) for the B n Fourier coefficients then becomes 


B 


n 


4 

7T 





A m sin m 6 sin n 6 d 6 


(A6) 


By interchanging the order of integration and summation, equation (A6) may be rewritten 

as 


B 


n 


I A m y ^ (x B - x^jsin mfl sin n8 dd 


m 


(A7) 


The type of polygonal wing planform initially to be considered is shown in figure 2. 

It is a tapered sweptback wing with straight leading and trailing edges and with stream- 
wise tips. In this case the geometric variables for complete planform definition may be 
expressed in three parameters: aspect ratio A, taper ratio X, and sweep angle A. An 
extension to more general planform shapes such as those encountered on a variable- sweep 
wing is given in the main text. 

For illustrative purposes, a field point will be considered where a Mach line through 
the field point intercepts both the leading and trailing edges. In this case, the field point 
is behind the trailing- edge root (x = c^j but forward of the leading- edge tip (x = kjy^ 
where c r = xq = kj, as shown in figure 2. The lower limit x^ in equation (A3) is given 
by the wing leading edge as kjy or -kj cos 9. The points where the Mach line x = xq 
through the field point intercepts the wing edges are designated as yj = -cos 6 * for the 
leading edge and y£ = -cos for the trailing edge. Consequently, for span stations 
inboard of 62 , the upper limit xg is c r - k2 cos 9, and for span stations between 
and 62 the upper limit is xq. After substitution of the appropriate limits, equation (A7) 
for the Fourier coefficients becomes 



sin m 6 sin n 6 d6» 


Upon integration and substitution of the relation 
xq = -kj cos 6 \ = c r - k 2 cos 82 

which is valid in the region c r g xq ^ kj, equation (A8) becomes 


kj cos 8 sin m# sin n 6 + xq sin mfi sin n6>)d0/ 


(A8) 


(A9) 
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APPENDIX A - Concluded 


Bn = ' i E A m k l< cos 6t f(0j) - g(*i) - cos % f(|)- G /tt 


m 


k 2< cos e 2 f (e 2 )- g(0 2 ) 


COS 0 2 Ffr- J - G 


(A10) 


where 


F(e) = | 


sin(n m )8 sin(n + m)fl 
n - m ” n + m 


(All) 


G(e) =\ 


sin(n - m + 1)0 sin(n - m - 1)0 _ sin(n + m + 1)0 _ sin(n + m - l)f 


n - m + 1 n-m-1 n+m + 1 


n + m - 1 


(A12) 

Simple indeterminate forms which are encountered in equations (All) and (A12) may be 
evaluated as follows : 


lim 
P ~ o 


sin p 6 
P 


= e 


(A13) 


It should be restated here that the B n Fourier coefficients obtained from equa- 
tion (A 10) are limited in application because of the assumptions inherent in the derivation. 
Equation (A10) is restricted to straight- edged sweptback wings with lateral symmetry and 
uniform chordwise load and is specifically formulated for field points in the region 
c r = xo = kj. These limitations are carried over to the evaluation of the C n 
coefficients in appendix B. 
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APPENDIX B 


EVALUATION OF THE C n COEFFICIENTS 

The C n coefficients are evaluated from equation (8) by using equations (A9) 
and (A10) and therefore are subject to the same restrictions as were the B n Fourier 
coefficients evaluated in appendix A. That is, they are for a polygonal planform (straight- 
edged) wing with lateral symmetry and for a uniform chordwise load. The field point is 
also located in the region c r = xq = kj. The extension to wings with less restricted plan- 
form and field-point location is given in the text. 

With the substitution of equation (A10) into equation (8), the expression for the C n 

coefficients becomes 



From equation (A9) it can be seen that 


dxg = k][ sin d0j \^ = 6 - (B2) 

and 


dxg = k£ sin #2 



Substituting equations (B2) and (B3) into the first and second integrals, respectively, of 
equation (Bl) gives the expression for C n 






(B4) 


Integration of equation (B4) yields 
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APPENDIX B - Continued 





wi> - 1 1 



where 



H(0) = - 


(n + m)(n + m + l)(n + m + 2) 

(n + m - l)(n + m)(n + m + 1) 

sin(n + m - 2)0 

2 sin(n - m)0 

(n + m - 2)(n + m - l)(n + m) 

(n-m - l)(n - m)(n - m + 1) 

sin(n - m + 2)0 

sin(n - m - 2)0 


(n - m)(n - m + l)(n - m + 2) (n - m - 2)(n - m - l)(n - m) 


(B5) 


(B6) 


1(0) = cos 0 


1 

i COS 

0 

- Gfli 

2 

\2/ 

\2J 


(B7) 


The simple indeterminate forms which are encountered in equation (B6) may be 
evaluated by using the form of equation (A13) as follows: 


lim 

(n-m) — 0 


sin(n - m)0 

(n - m - l)(n - m)(n - m + 1) 


e 


lim 

(n+m-2) - 0 


sjji(n + m - 2)6 

(n + m - 2)(n + m - l)(n + m[_ 


6 

2 


lim 

(n-m +2) — 0 


sin(n - m + 2)6 

(n - m)(n - m + l)(n - m + 2)_ 


e 

2 


lim 

(n-m-2) 


- 0 


sin(n - m - 2)0 

(n-m - 2)(n - m - l)(n - m) 


e 

2 


lim 

(n-m) - 0 . 


sin(n - m + 2)6 


sin(n - m - 2)6 


(n - m)(n - m + l)(n - m + 2) (n-m - 2)(n - m - l)(n - m) 


= 0 cos 2 9 - sin 2 0 
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APPENDIX B - Concluded 


Once again it should be stated that the limitations inherent in the B n Fourier coef- 
ficients represented by equation (A10) also apply to the C n coefficients represented by 
equation (B5). The derivation of these equations , however, is indicative of the methods 
used for extension to less restricted planform and field-point locations. The limitation 
to uniform chordwise load, however, is retained throughout the analysis. 
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APPENDIX C 


ANALYTICAL SUMMATION OF SERIES 

The series summations over n of equations (6) and (7) with the use of equa- 
tions (A 10) and (B5), respectively, are not rapidly convergent and generally require a 
large number of terms to be summed before a satisfactory numerical result is achieved. 
With equation (B5) for the C n coefficients substituted into equation (7) and with the order 
of the summations over n and m interchanged, the summation over n may be 
obtained independently of the summation over m. An example of the methods used to 
obtain one of the summation terms is presented in appendix D. Similar methods are 
applicable to the other summation terms. The results of the summation are summarized 
as follows: 


n sin n^jcos 0 F(0) - G(<9)] = D m (0) 


n=l 

CO 


Y, n sin n#o F (f) = "^m^) 
11=1 

oo 

X nsinn0 o = -E>m(f 

n=l 

oo 

Y n sin n% H(0) = E m (0) 


(Cl) 


(C2) 


(C4) 


where 


D m (e) = y cos e {\ cos m6 o lQ g 


sin 2 (0o + 0) 


m+1 

2 


r=l 


sin^(0Q " 0) 

(m cos m0Q cos Oq - sin vciOq sin log 


A . n n 2 sin sin(m + 1)0 

2 V sm 2r 6 sin(m : 2r \0 + _ Q ^ + 0 sin m0 Q ( 

Lt 2r m + 1 

sin(2r - 1)0q sin(2r - 1)Q\ 


m+1 
2 2 




2m 


isin m 6 sin 6q - (m sin m0Q cos Oq + cos m0Q sin 


sin(0Q + 0) 1 - cos(0q - 0) 

[sin(0Q _ e ) 1 “ cos(0q + 0)] 

m+1 

cos(2r - 1)0q sin(2r - 1)0 


2r - 1 


r=l 


2r - 1 


(o t %) (C5) 
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APPENDIX C — Continued 


;%W = ~ sin m0 


K 0) - g m(cos 2 0 + cos 2 cos m0Q - 2 cos 6q sin 0q sin m0Q \ — log — — Z® 

L X ' J 2 [sin 2 ( 6q - 6) 


Y sin 2r0 o sin 2r0 \ j 

4 Z o 1 + 4 cos e K n cos m0 O cos e 0 - sin m0 o sin 0 Q ' 


i> , fi - c ° s («o - «> 


x\±log ^ +log- x u 7 +4 

) 2 _sin 2 (0 o -6) g 1 - cos(0q + 0) 


sin(2r - 1)0q sin(2r - 1)0 


+ | m sin m0 o (cos 2 0 + cos 2 0q) + 2 sin e Q cos 


x \26 + 4 


cos 2r&Q sin 2r 
2r 


\ cos 0(m sin m0Q cos 6 »q + cos m0Q sin 


m+1 
~ *r~ 

A Y cos(2r - 1)0q sin(2r - 1)0 

Z 2r - 1 

r=l 


sin(m -1)0 sin 0g sin(m + 1)0 sin 30^1 j 


olm - 2) + — 

o m - 1 m + 1 


2 m 


■cos 0 sin m0Q sin 


m{2 + cos 20) . 


| lsin(m + 1)0 sin 0q 


(# £ 0q) (C7) 



APPENDIX C - Concluded 


E m (0) = j cos 0o(m cos 0Q cos m&o - sin 0q sin log (cos 2 0q) " 4 ^ 


m+1 . . •*. „ 9 a 

V (-1) sin r 0 q 


+ ^ cos 0 q(hi cos sin m0Q + sin 0 q cos 20q + 4 


1 sin(m - l)0o sin 0 q sin(m + 1)0 q sin 3 0 q 

+ - (m - 2) + 

8 m - 1 m + 1 


m+1 . „ \ 

V (-1) 1 cos r0Q sin r0Q \ 


(m - m(2 + cos 20 q) 

iS — ±1 cos 0n sin m0n sin 0n + sin(m + 1)00 sin 6q 

2m u u j 4(m + 1) 


(0 = 0Q) 


and the indeterminate form in equations (Cl) and (C8) may be evaluated as follows : 


sin(m - 1)0 O sin 0 O . Q 

lim = 0o sin 0 q 

m — 1 m - 1 J 


D m (0) E m (0) 

To assist in numerical calculations, the functions : — and : — are 

m sin m0Q m sm m0Q 

presented in graphical form in figures 8 and 9 for odd values of m up to m = 39. A 
curve for m approaching 00 has been added to show the limit toward which these func- 
tions tend as m grows very large. Figures 8 and 9 have been prepared for- fixed values 
of yo corresponding to values of 6q at 7.5° increments from 90° to 172.5°. The 
curves at 0q = 120° for the odd values of m which are multiples of 3 have been deleted 
from figures 8(e) and 9(e). They are presented separately in figure 10 for D m (0) and in 
figure 11 for E m (0). These extra figures are necessary because sin m0o equals zero 
in the denominator of the functions presented in figures 8 and 9 at 0o = 120° for these 
values of m. 
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APPENDIX D 


DERIVATION OF ANALYTICAL SUMMATION OF SERIES 

The derivation of the closed- form result (eq. (C5)) of the summation expressed in 
equation (Cl) is shown in detail. Equation (Cl) is rewritten 

OO 

D m (0) = ^ n sin n6 0 [cos 0 F(0) - G(0^ (Dl) 

n-1 


The substitution of the expression for F(0) from equation (All) and for G (0) from 
equation (A12) into equation (Dl) yields 




1. sin(n - m + 1 
4 n - m + 1 


sin(n - m - 1)0 sin(n + m + 1)0 


m - 1 


n + m + 1 


sm(n + m - 1 
n + m - 1 


Equation (D2) may be simplified to 



As previously mentioned, if the wing has lateral symmetry, then only odd values of n 
occur. Consequently, with the substitution of n = 2r - 1 where r takes on consecutive 
integer values, equation (D3) may be rewritten as 

oo 

D m (0) = Yj sin ( 2r " !) e 0 

r=l 


si n(2r - xn - 1)0 sin(2r + m - 1)6 
2r - m - 1 2r + m - 1 
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APPENDIX D - Continued 


After rearranging and combining terms, equation (D4) becomes 


D m ( 0 ) = Tf cos < 9 <cos m00 


cos 2 t{ 0 q - 0 ) 


r=l 


2r 


cos 2t(6q + 6) 


m+1 

T 


r=l 


2r 


- y 


r=l 


cos 2r(0Q - 6 ) 


2r 


m+1 

2 

Y cos 2r(e 0 + 6 ) cos(m + 1 ){ 6 0 - 6 ) cos(m + 1 )( 6 q + 6 ) 


r=l 


2r 


m + 1 


m + 1 


+ sin m0Q 


m+1 m+1 

___ — 

sin 2r(% + 6 ) ^ sin 2r(0o - 0 


r=l 


2r 


r-1 


2r 


sin(m + 1)(0 q + 0 ) sin(m + l)(0p 
m+1 m+1 



cos ui 6 q cos 6 q - sin hi< 9 q sin Oq 


n 


r=l 


cos(2r - l)(^o - 0 ) 
2r - 1 


m+1 

OO -J- 

1 V cos(2r * l)( 6 >o + 0) V sin(2r - 1)00 sin(2r - 1) 

2 Z 2r - 1 l-j 

r=l 


r=l 


^m sin m0o cos 6 q + cos m0Q sin 00^ 


2r - 1 
m+1 


- n ~ sin m0 sin 0 n , 

2m u 


cos(2r - 1)00 sin(2r 


r=l 


2r - 1 


(D5) 


The following summation Is presented In reference 13: 

1 


cos r6 1 , 

“T— = 2 log 


r=l 


2(1 - cos 0) 


(0 < 0 < 2 it) (D6) 


From equation (D6), the following summation may be obtained : 

1 


r=l 


cos 2r 0 1, 

— 2F- = 4 l0g 


2(1 - cos 20) 


(0 < 0 < it) (D7) 
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APPENDIX D — Continued 


Equation (D7) may be rewritten as 


S2|i£« = - I log(4 sin2<J) 

r=l 


(0 < 6 < n) 


The series of equation (D6) may be written as the sum of two series 


(D8) 


cos rd _ V cos(2r - 1) 6> y 

r L 2r - 1 + L 


cos 2r 6 
2r 


r=l r=l r=l 

The substitution of equations (D6) and (D8) into equation (D9) yields 

oo 

y cos(2r - 1)0 _ 1 i / sin 8 x 2 

A 2r - 1 4 ° \l - cos 8, 

r=l 


(0 < 6 < 77 ) 


(D9) 


(DIO) 


The cosine series of equations (D8) and (DIO) and the functions they represent are 
even functions of 6 and are periodic in intervals of 2 it. Consequently, the intervals of 
convergence of equations (D8) and (DIO) may be written 


-7T < 8 < 0, 0 < 6 < it, 7T < 8 < 2ir 


It is therefore seen that, in these two equations, the series converges to the function 
everywhere except at 6 = 0 and at 6 = ipir where p is any integer. 

With the use of equations (D8) and (DIO), the infinite series in equation (D5) may be 

written as 


l 


cos 2r(6> 0 - 6) 
2r 


l 

T* ~ 


cos 2r(#Q + 6) 

__ 



sin 2 (0Q + 6) 
sin 2 ( 6q - 6)_ 


(Dll) 


r=l 


cos(2r - 1)(0 q - 6) 
2r - 1 


-l 


r=l 


cos(2r - l)(6>o + 6) 
2r - 1 


3 10 * 


sin( 6q -6) 1 - cos(0o + 6) 


sin(0o + 6) 1 - cos{6q - 6) 


(D12) 


In equations (Dll) and (D12), the range of the angle arguments are (since, in general, 

^ ~ it) - ■§ = (6q - 6) =?- and ir = (6q + 6) £ 2tt. Therefore, the use of equations (Dll) 

and (D12) in equation (D5) is valid everywhere except at 6 = 6q. 
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APPENDIX D - Concluded 


Substitution of equations (Dll) and (D12) into equation (D5) yields, after 
simplification. 


m+l 


D. 


(B) =y cos 6 < 


m 


cos m% log 


sin^(^o + 6 ) 


sin^(0Q - 6 ) 


2 V sin 2r 0 sin(m ~ 2r )6 + 2 sin % sinym + i, ; ^ 

^ 2r m + l 

r=l 


: sin m t 


+ ^m cos xxiOq cos - sin m 0 Q sin 0 qK log 


sin(% + 6 ) 1 - cos(0q ” 0) 


sin(0Q - 0 ) 1 - cos( 0 o + 0 ) 


m+l 


r=l 


sin( 2 r - 1)0q sin( 2 r - 1 ) 6 \ 
^ 2r - 1 


sin m 6 sin 
2m u 



sin m Oq cos 0 q + cos m ^0 s ^ n 


%) 


m+l 



cos( 2 r - 1 )< 9 q sin( 2 r - 1)0 
2r - 1 ~~ 


(6 8 q ) (D13) 


Equation (D13) is the same as equation (C5). The other summation expressions which are 
represented by equations (C7) and (C8) may be obtained in a similar manner. 
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APPENDIX E 


CALCULATION OF THE A m COEFFICIENTS 


The desired wing span wise load distribution is considered to be known so that, for a 
given planform, the spanwise distribution of section lift coefficient c z is also known. If 
an analytical expression for the c^ distribution is known, then the A m coefficients of 
equation (11) can be found by a Fourier expansion. However, if the c^ distribution is 
known only numerically, then the A m coefficients may be approximated by simple 
numerical methods for all but the most irregularly shaped c i distributions. Two such 
methods are the 17-point odd harmonic numerical method for the sine terms presented in 
reference 14, and the method in reference 11. 

The A m coefficients for elliptical spanwise load (cc^) distribution on straight- 
tapered wings with streamwise tips may be obtained as follows: 


A 


m 


1 C 71 4 f 71 / 2 

= — \ c; sin m0d0 = — \ c ? sin m6 d& 

vJ-v 71 J 0 1 


where, for elliptical spanwise load 


(El) 


cc l = i c r(l + *)C L sin 6 


(E2) 


In the range of 6 from 0 to the chord distribution may be expressed as 

u 


c = c r jl - (1 - X)cos o] 

Substituting equations (E2) and (E3) into equation (El) yields 

'7r/ 2 
'0 


o ^ e s 


A m =i(l + A)C L £ -Si 

71 z -'o 1 - 


sin 6 sin m 0 


dd 


(1 - X)cos 8 

The solution for the integral of equation (E4) can be written in the recursion form 


(E3) 


(E4) 


L m 4 1 + X 


c L n-x 


^m-l “ ^m+1 


(E5) 


and the recurring K m functions are evaluated 
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APPENDIX E - Concluded 


K 0 = 


X 


rf (2 - X) 


f i r . . 

2 + sm 


K 


1 - 


K 0 
1 - X 


1 

2 


XT , jr _ • H17T ^ it 

K m _ 1 + K m+ 1 - - — sin 2 + J x K m 

Table 1 presents values of A m /c L for values of X from 0 to 1.0. The values 
of A m /CL for X = 0 are for a wing with pointed tips and do not represent a realis- 
tically possible lift coefficient distribution. They are included in the table only to indi- 
cate limiting values of A m /C L as X - 0. In general, the first few entries in the 
table (at the lowest values of m) were obtained with the recursion formula of equa- 
tion (E5). At the higher values of m, the table entries were obtained by numerical 
methods. 

The method of solution for camber presented in this report is not well suited for 
use where a uniform spanwise section lift coefficient distribution is desired. However, 
approximate results can be obtained with the A m j Cl values presented in table 2 which 
were computed by the method of reference 11. The values of A/nm obtained by 
Fourier expansion for a uniform c i distribution are included in the table for comparison 
purposes. 
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TABLE 1,- A m /C r FOR ELLIPTICAL SPANWISE LOAD DISTRIBUTION ON STRAIGHT -TAPERED WINGS 
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Figure 1.- Sketch of wing showing the partial wing loading on the region S*. 
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Figure 7.- Comparison of experimental and calculated twist distribution for 
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